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Abstract. As a preparation for the dynamical investigations, this paper begins with 
a short review of the three-layer gravastar model with distinguished attention to the 
structure of the pertinent parameter space of gravastars in equilibrium. Then the 
radial stability of these types of gravastars is studied by determining their response for 
the totally inelastic collision of their surface layer with a dust shell. It is assumed that 
the dominant energy condition holds and the speed of sound does not exceed that of 
the light in the matter of the surface layer. While in the analytic setup the equation of 
state is kept to be generic, in the numerical investigations three functionally distinct 
classes of equations of states are applied. In the corresponding particular cases the 
maximal mass of the dust shell that may fall onto a gravastar without converting 
it into a black hole is determined. For those configurations which remain stable the 
excursion of their radius is assigned. It is found that even the most compact gravastars 
cannot get beyond the lower limit of the size of conventional stars, provided that the 
dominant energy condition holds in both cases. It is also shown — independent of 
any assumption concerning the matter interbridging the internal de Sitter and the 
external Schwarzschild regions — that the better is a gravastar in mimicking a black 
hole the easier is to get the system formed by a dust shell and the gravastar beyond 
the event horizon of the composite system. In addition, a generic description of the 
totally inelastic collision of spherical shells in spherically symmetric spacetimes is also 
provided in the appendix. 

PACS numbers: 04.40.Dg, 97.60.Lf, 04.20. Jb, 95.36.+X 
1. Introduction 



There are more and more astrophysical observations justifying the existence of extremely 
compact massive objects with size close to their Schwarzschild radius [H EJ |3]. It 
is widely accepted that these observations also provide indirect justifications of the 
existence of black holes (BHs). Nevertheless, there are also alternative ideas about 
claiming that exotic states of matter may exist which could stabilize extremely compact 
stars suiting to the aforementioned astrophysical observations (see, e.g., [4] for a recent 
review). One of the most popular among these types of BH mimicking objects is the 
gravitational vacuum star (gravastar) model which has received considerable attention 
not least because its relation to the concept of dark energy. In this model of Mazur and 
Mottola [5] an interior de Sitter spacetime region is connected via three intermediate 
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layers to an outer Schwarzschild solution such that the radius of the outermost layer is 
supposed to be slightly larger than the Schwarzschild radius of the system. 

It is worth mentioning that in advance to the gravastar model there were several 
constructions in which the matching of a de Sitter region to the Schwarzschild spacetime 
was applied. For instance, to get rid of the r = singularity of the Schwarzschild 
spacetime, Frolov, Markov and Mukhanov in [6] proposed a matching of a de Sitter 
interior to it at a small radius of the Planck scale ensuring thereby that the curvature 
remains bounded everywhere in the yielded spacetime. Dynamical investigation of this 
model was already carried out in [7]. 

Note that up to certain extent the model used in [6] could be considered as the 
precursor of the gravastar model — although in the latter the matching was made in 
a more elaborated way — and the outermost matching surface, at the boundary of the 
Schwarzschild region, was supposed to be arranged such that its radius is slightly larger 
than the pertinent Schwarzschild radius. 

Once such a model is set up the following questions manifest themselves: 

(1) What type of physical process may produce such a gravastar? 

(2) Is a gravastar stable? 

(3) If it is, does it provide a viable alternative to BHs? 

While the first question has not been tackled yet even the second question turned 
out to be too complex within the original model of Mazur and Mottola — although in [5] 
an argument claiming for its thermodynamical stability was given — as it is composed 
by making use of three different types of regions with unspecified matter. To reduce the 
related ambiguities Visser and Wiltshire [8] introduced a simplified three-layer gravastar 
model where the interior de Sitter region is matched to the exterior Schwarzschild 
spacetime via a single matter shell. This model is simple enough to carry out various 
analytic investigations by making use of the thin-shell formalism of Israel j9]. Visser 
and Wiltshire besides deriving the basic relations determining the evolution of gravastars 
also carried out the first investigation of their radial stability. Since then the stability 
of gravastars has been studied by several authors within this simplified model or within 
its continuum correspondence [10] . Results relevant for radial stability can be found in 
[TT|[T2"]. and in case of electrically charged gravastars in [13]. The stability has also been 
investigated with respect to axial perturbations [IH fT5] . 

In all of these investigations attention was restricted to the space of gravastars in 
equilibrium, i.e. the radial stability was investigated by determining the response of a 
gravastar to a slight formal change of the underlying effective potential. For instance, in 
[T6| [T7] the excursions of gravastars was investigated in such a way that their evolution 
started with carefully prepared initial conditions. In all of the pertinent investigations it 
was demonstrated that by suitably adjusting the equation of state (EOS) of the matter 
forming the surface of the gravastar the subspace of 'stable' gravastars may always be 
ensured to be of non-zero measure. Nevertheless, it was also found that whenever the 
measure of the subspaces of the configuration space representing stable and unstable 
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gravastars is compared the former is always found to be negligible with respect to the 
latter. This observation was commonly interpreted that gravastars may not offer a 
viable alternative to BHs. 

The main purpose of the present paper is to determine the response of a gravastar 
in equilibrium to the arrival of a dust shell onto its surface. This is done not merely 
by considering some formal change of the effective potential determining the state of a 
gravastar but also by making use of the full dynamical setup. For the sake of definiteness, 
we assume that the surface of the gravastar and the dust shell collide in a totally inelastic 
manner. In addition, concrete EOSs are chosen and it is assumed that the dominant 
energy condition (DEC) holdfl and the speed of sound in the surface of the gravastar 
does not exceed that of the light. Then the relevant non-linear problems — the basic 
equations of which are based on the dynamics of spherical shells — are solved by using 
analytical and numerical approaches. In this way we can not only study the excursion 
of particular gravastar models but also determine the maximal mass of the dust shell 
colliding with the surface of the gravastar without converting the latter into a BH. 

This paper is organized as follows. In section [2] some of the basics of the Visser 
and Wiltshire three-layer dynamics gravastar model are recalled using dimensionless 
variables. As a preparation for the aforementioned dynamical investigations a short 
survey of the configuration space of stable gravastars is also provided. (Although there 
are no completely new results in this section, we believe that this review provides a 
good reference frame for the results of the succeeding dynamical investigations.) In 
Section [3] the dynamics of the system composed by the spherically symmetric dust shell 
falling onto a stable gravastar, along with their collision, is described. Section @] reports 
about our analytical and numerical results concerning the dynamics of maximally loaded 
gravastars, while Section [5] contains our concluding remarks. Finally, in the appendix, 
a generic description of the totally inelastic collision of spherical shells in spherically 
symmetric spacetimes is provided. Throughout this paper the geometrized units, with 
G = c = 1, are applied. 

2. Gravastar model of Visser and Wiltshire 

Throughout this paper considerations will be restricted to the three-layer spherically 
symmetric gravastar model of Visser and Wiltshire [8]. This simplified model consists 
of an external Schwarzschild vacuum region with the mass parameter M, representing 
the total ADM mass of the gravastar, an interior de Sitter region with energy density 
Po, and an infinitesimally thin shell at the mutual boundary of the aforementioned 
two regions at radius a with surface energy density o and surface tension 9. The line 
element representing the interior de Sitter and exterior Schwarzschild spacetimes is given 
in the form of (lA.lj) and (1A.2I) . where the pertinent forms of Mk(t) are specified in the 

1 Recall that the dominant energy condition guarantees that the concept of causality is properly 
adopted in general relativity. Therefore, the exclusion of matter models yielding its violation seems to 
be preferable. 
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paragraph below (1A.2D . Note that the radius is allowed to vary in time. As the shell is 
assumed to coincide, at any moment, with one of the SO (3) group orbits, the quantities 
characterizing the shell are only the functions of the radius. To have a closed system 
of equations governing the evolution of the shell we also need to assume the existence 
of an EOS. As was noted above a and 9 must be functions of the radius exclusively. 
Thereby, in virtue of the implicit function theorem, the EOS, relating the tension, 
9, of the surface to the mass density, a, of the surface of the gravastar, may always be 
assumed to possess the form 9 = 9(a). This section, besides recalling the most important 
equations determining the dynamics, also provides a short comprehensive review of the 
basic properties of the above-outlined three-layer gravastar model. 



2.1. Parametrization 

In studying the gravastar model the use of the following dimensionless parameters turned 
out to be advantageous. They are the dimensionless radius a = a/(2M), surface density 
S = o~M, tension O = 9M and the dimensionless de Sitter density parameter, r\ = 8kM 2 , 
where k = 47rp /3§| Note that in the chosen setup positive tangential pressure does 
correspond to negative surface tension. To avoid the appearance of the event horizon 
in the Schwarzschild region the inequality 1 < a is required to be satisfied. Similarly, 
in order to exclude the presence of a 'cosmological' horizon in the de Sitter region the 
inequality a < rf 1 ! 2 is imposed. Accordingly, in case of the considered gravastar models 
the value of r\ is restricted to the interval < rj < 1. Finally note that in dynamical 
situations the time dependence may conveniently be expressed via the dimensionless 
proper time measured by observers moving radially along the world-sheet of the shell, 
in 2M units. 



2.2. Dynamical equations 

The basic equations determining the dynamics of the shell of the gravastar can be derived 
in various ways. The most widely applied method uses the thin-shell formalism of Israel 
[9] based on the results given by Sen, Lanczos, Darmois [HI EE] [20] . In this approach the 
induced metrics, h AB and h AB , on the mutual boundary of the two spacetime regions— 
the inner and the outer — on the two sides of the shell are assumed to coincide, i.e. the 
metric on the shell can be given as Uab = h AB = h AB , while the discontinuity of the 
pertinent extrinsic curvature tensors, K\ B and K AB , are related to the surface-energy 
density Sab = diag(a, —9, —9) via the matching conditional 

K\ B - K AB = -8tt {S AB - h AB (h CD S CD )} ■ (2.1) 

2 Whenever r\ = k = the systems degenerates to a shell moving in a Minkowski-Schwarzschild 
spacetime. 

3 This relation corresponds essentially to the Lanczos equation. 
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By making use of the pertinent results of [8], along with the above-introduced 
parameterization, these junction conditions can be seen to take the form 




1 - r]a 2 + a 2 - Vl-o^ + d 2 , (2.2 



. 1 — 2na 2 + a 2 + aa I — a 1 /2 + a 2 + aa\ 

~ } = ^~ ! Vi-^ + d* vi-^ + d 2 ~ J ' ( } 

where the over-dot denotes the derivative with respect to the dimensionless proper time 
along the shell. 

As the energy density of the shell is always assumed to be non-negative, the first 
square root on the left-hand side of (12. 2p must be greater than the second one, which 
implies that a < r? -1 / 3 . Note that this restriction, viewed as a condition on the value of 
a, is more demanding than the previous one excluding the presence of a cosmo logical 
horizon in the inner de Sitter region. In virtue of (12. 2p it is straightforward to see that 
whenever E is non-negative it attains its maximum value for vanishing a for any fixed 
value of a. 

In describing the dynamics of gravastars it turned out to be advantageous to 
introduce an 'effective potential' defined via the relation 

a 2 

V{a-r 1 ) = -—. (2.4) 

Assuming that E is positive, by a straightforward algebraic manipulation, one finds that 
V can be given as a function of a and E = S(a) — it also depends on the dimensionless 
de Sitter density parameter rj — and it reads 

2" 



V(a,Z(a);r)) = - 



G4ir 2 a E 2 (a) V 167ra 2 E(a) 



, (2.5) 



which, in virtue of (12. 4 p — despite its complicated analytic form — has to be smaller than 
or equal to zero. 

Once E is known as a function of the dimensionless radius, a, so is V(a;r]), and, 
in turn, (12 .4p can be used to determine the motion of the shell. The next subsection 
shows that the EOS of the shell, along with the conservation law, D Sab = — where 
Da denotes the covariant derivative operator compatible with the metric, hAB, on the 
shell — can be used to determine the functional form of E(a). 



2.3. Equation of state of the shell 

Recall first that the EOS of the surface matter of the gravastar, if given in terms of the 
dimensionless variables, may be assumed to possess the form O = ©(E). Admittedly, 
the entire concept of an infinitesimally thin shell may be considered to be somewhat 
artificial. Nevertheless, in our investigations we need to apply some assumptions on the 
type of its matter. The use of a suitable energy condition is essential. In this respect 
the use of the DEC seems to be the most appropriate which is supposed to guarantee 
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that the concept of causality is properly adopted in describing the evolution of coupled 
gravity matter systems. Note that the less stringent weak and null energy conditions are 
automatically satisfied whenever the DEC is guaranteed to hold. For an infinitesimally 
thin shell the DEC can be seen to hold whenever |0| < E [21] . 

Throughout this paper the function O = ©(E) will be assumed to be continuous 
and piecewise differentiable. To avoid dynamical instabilities we shall also assume that 
the square of the speed of sound, Cg = — d©/dE, is non-negative, and to be compatible 
with the concept of relativity that cj is required to be less than or equal to the square 
of the speed of light. Accordingly, unless otherwise stated, it will be assumed that 
< c 2 s < 1. 

It is well-known that the radial conservation equation is 

i M = e i M • < 2 - 6) 

where r denotes the dimensionless proper time along the radial trajectories generating 
the shell. (12 .6p may be derived either by referring the conservation law D a Sab = or 
by making use of equations (12.21) and (12.31) as in [8]. By introducing the dimensionless 
radius a, instead of r, as an independent variable (12. 6p takes the form 

a^- = 2(9 - E) . (2.7) 
da 

This is the very point where the EOS of the shell comes into play since by substituting 
the relation O = 0(E) into (12.71) the functional form of E = S(a) can be determined. 
Note that as 9 - E < 0, in virtue of ( 12. 71) . S(a) has to be a monotonic decreasing 
function. Since (12. 7\\ is a separable differential equation, an implicit solution of it can 
be given as 

- expfir^, (2.8) 



a \2Jz 9(E) 

where the integration constant So is fixed such that So = E(«o)E 

To have some specific examples in later sections we shall use the following 
three types of EOS. The simplest EOS — that has been used in most of the former 
investigations of gravastars — is linear homogeneous possessing the functional form 
©(E) = — wq E, which, however, is not satisfied for a generic linear functional relation. 
In this case cl = wo and 

/ a \ -2(i+too) 

E(a) = E - . (2.9) 



Clearly, it is of interest to consider more complex systems, as well. Since dO/dE 
is related to the speed of sound and its value at an equilibrium point will play an 

4 Throughout this subsection the integration constant are determined by referring to certain ao value. 
Although latter ao will label equilibrium states, it should keep in mind that «o does not have such an 
interpretation yet. 
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important role in the stability of the gravastar, it seems to be reasonable consider the 
EOS with varying sound speed at the equilibrium point. The simplest generalization 
of a homogeneous linear EOS would be the non-homogeneous linear one. However, we 
would also like to have the DEC to hold which requires O to vanish at E = 0. Therefore 
we use instead the 'broken linear' EOS 

Here the parameter Si signifies the value of the breaking point in the E range, while for 
the values of the slopes, w± and w 2 , before and after the matching the inequalities 
< u>i,u>2 < 1 hold. Then ( 12.8f) can be evaluated explicitly. In doing so the 
dimensionless radius of the matching point «i, gets to be determined by the relation 

E x = as 



^l/(2»i;i+2) 



i «b(g) ,if£o<£i 

ai ~^ ^o(i+ W2 )+s 1 ( Wl - W2 )A 1 /( 2 ^+ 2 ) - U ' U) 



Then, whenever E < Ei 



w2-wi+(i+wi)(^y 



2tu 2 +2 



E fa) , if a > ai , 



while for E > Ex 



s<..)= ; . . te+J — — ' ,fa - ai .2.13) 

Si f^) , if a > Ul 



Another obvious choice is a polytrop EOS. Note, however, that in general a polytrop 
EOS is supposed to relate the rest energy density e, as opposed to the total energy 
density a, to the pressure in which case it takes the functional form p = Ae K , with 
parameters A and k. Thereby, whenever we replace e by a the 'inverse form' of the 
polytrop EOS reads as (see, e.g., [22]) 

/ 0\ 1/ft O 



A) k-1 
The solution of (12. 8p can be given explicitly as 



where O is related to the value of E by (12.141) . 

It is straightforward to check that the DEC is satisfied if and only if k G (1,2]. 
Assume, now that if we are given both of the values Eo and ©o and the EOS is chosen 
to possess the form (12. 14ft then k has to take value in the interval [1 + Wq,2\, where 
w = — ©o/S . Then the square of initial speed of sound at E can be given as 
Cslso = kw /(1 + w ). 



-7- 



Merse Elod Gdspdr and Istvdn Rdcz 

In closing this subsection let us emphasize again that whenever an EOS, G = 0(E), 
is chosen E may be given as a function £ = E(a) and, in virtue of ( 12. 5p . the effective 
potential can be given as an explicit function V(a; rf) depending exclusively on the 
dimensionless radius a and the dimensionless de Sitter density parameter rj. 



2-4- Gravastars in equilibrium 

In this subsection considerations will be restricted to the case of gravastars in 
equilibrium. Equilibrium, as usual, signifies those configurations in rest for which forces 
are also balanced, i.e. V(ao; rf) — and V'(ao; rf) = 0, where «o stands for the value of 
a in equilibrium, while prime, denotes derivative with respect to a. As is expected the 
vanishing of V and V is equivalent to the vanishing of a and a, which can be verified by 
making use of (12.41) . In depicting of the configuration space of gravastars in equilibrium 
we shall apply the values E = E(a ) an d O = Q(a ), which, by making use of (12. 2|) 
and (12.31) . along with the vanishing of a and a at « = «o, can be given as 

° 87TQ!o 

and 



1- V a 2 - VI -«o (2-16) 



8 = ^- _ I- w y (2 . 17) 

Using the algebraic equations of ( 12.161) and ( 12.171) and the condition that Eq(?7, cxo) > 0, 
it can be shown that for a gravastar in equilibrium the surface tension, ©0(77, «o), is 
always negative. This then justifies that no exotic matter (with negative pressure) is 
required to form the shell of a gravastar in equilibrium. 

Note that whenever an EOS is given only one of the values of rj and can be 
chosen freely because, in virtue of ( 12. 16ft and 1 12. 17ft . w$ = — Oo/S has to be consistent 
with the EOS. On the other hand whenever no definite choice is made for an EOS then 
the configuration space of gravastars in equilibrium may be parameterized by rj and a^. 
In this case the value of w (r], a ) gets to be determined by relations (12.161) and (J27TTJ) . 

In depicting the configuration space of gravastars in equilibrium it is useful to 
identify the region on the (ceo,r)) plain where the DEC holds (see figure [2]). In doing so 
recall first that whenever E is non-negative for a gravastar in equilibrium O has to be 
negative. This immediately implies that the DEC reduces to < Wq < 1; furthermore 
the upper boundary of the subregion of the (cto,^) plain where the DEC is satisfied 
is signified by the wq = 1 curve. This curve — determined by the implicit relation 
w (77, a ) — 1 — can be given as 



DEC 60a§ - 36ct - 25 + (5 - 6a ) \/l00a% - 124a + 25 
W M = 128 08(00-1) ' (2 - 18) 

which intersects the a -axis at the value a = 25/24 (see figure [2]). Consequently, 
whenever the DEC holds «o has to be greater than or equal to 25/24. This value of the 
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minimal radius immediately imposes a strict limitation on the compactness, expressed 
by the 2m jr ratio, of a gravastar in equilibrium. It also follows from d 2 . 1 8 [) that for 
small 7] the approximation r)^£ ~ |«o 3 holds. 

For comparison it is useful to recall that the value of the ratio 2m/ r is also restricted 
in the case of conventional star models composed by matter satisfying the DEC. In the 
case of spherically symmetric static configurations under various assumptions concerning 
the energy density and pressure of the star the ratio 2m jr has been investigated. A 
comprehensive summary of the pertinent results can be found in [23]. It is shown there 
that if pressure is non-negative (although it may be anisotropic) and the DEC holds, 
then the numerical value of the sharp upper bound for 2m jr is 0.963@ It is interesting 
that the corresponding upper bound for a gravastar in equilibrium is a little bit smaller 
24/25 = 0.96. One of the readings of this finding is that whenever the DEC holds the 
radius of the most compressed gravastar in equilibrium is definitely larger than the radius 
of the most compressed conventional star. This indicates then that instead of asking 
whether these gravastars can be distinguished from BHs, by any sort of astrophysical 
observations, it is more appropriate to ask whether they can be distinguished from the 
most compact stars. 

2.5. Stable equilibrium 

In deciding whether a gravastar is in equilibrium there has been no need to choose 
a specific EOS. Nevertheless, whenever we want to study the radial stability of 
these configurations we shall need to know the EOS, at least, in a sufficiently small 
neighborhood of So as merely the knowledge of ©o and So will not suffice then. 

Before turning to the radial stability of gravastars there are some conceptual issues 
to be mentioned. In spite of the familiarity of the dynamical equation a 2 /2 + V(a; rf) = 
it is better not applying automatically the analogous one-dimensional classical 
mechanical arguments. The main reason is related to the fact that V(a; if) is not an 
external potential. In the case of the classical mechanical problem the 'total energy' 
need not to be zero so we may simply either increase the kinetic energy term or change 
the potential energy by moving out of our system from its equilibrium. In both cases, 
since the potential is supposed to be intact, the total energy is changed. As opposed to 
this in the case of a gravastar if the initial velocity is changed the potential has also to be 
changed since the formal 'total energy' needs to remain zero. The change in the initial 
velocity a affects the initial surface energy density, S , therefore the solution S(a) of 
( 12. 8ft . and, in turn, V(a, T,(a);rj) is also affected. This verifies that it is better to avoid 
thinking of V(a, S(a); 77) as an external potential. 

In carrying out the radial stability of gravastars we shall assume that V smoothly 
depends on its indicated variables. This implies then that under a sufficiently small 
perturbation of a gravastar in equilibrium the sign of the second derivative, V"(ao), will 
be intact. In particular, since the perturbation increases the kinetic term, the potential 

5 This numerical value is claimed to be precise up to three digits in |23) . 
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has to become negative in a sufficiently small neighborhood of ao whence the yielded 
motion of the gravastar will be bounded if the value of V"(a ) is guaranteed to be 
positive. Accordingly a gravastar in equilibrium, with radius ao, is considered to be 
stable if V"(a ; 77) > 0. 

As was already mentioned, the vanishing of V(ao) and V'(ao) — in the classical 
mechanical analogy — corresponds to vanishing speed, a = 0, and, in turn, the forces 
are balanced. The equilibrium configurations can be classified according to the convex- 
concave character of the potential, i.e. it is reasonable to distinguish the following three 
categories. 

(a) Stable: A gravastar in equilibrium is considered to be stable whenever small 
perturbations of the system yield small local oscillations around the equilibrium 
state. In this case the potential has to be a convex function near ao. This may 
happen whenever V"(a ) > 0. (See, e.g., the transition from potential 1 to potential 
2 in figure [T]) . 

(b) Locally unstable: This occurs if the potential is concave (at least from one side) 
near ao. This means that small perturbations may yield non-local motions. This 
locally unstable situation may arise, e.g, whenever V"(ao) < 0. (For an illustration 
see graph 3 in figure [U where a± takes the role of ao). If V" is not continuous at 
ao but has well-defined left- and right-hand-sided limit values the equilibrium state 
will be called locally unstable if either of these limit values is negative. 

The vanishing of V"(ao) may also occur. In this case the sign of the least non- 
zero a-derivative of V determines the character of the potential. For instance, the 
graph of V may change its character from concave into convex or vice versa through 
an inflexion point whenever V" is the first non-zero higher order derivative. This 
means that the equilibrium is locally unstable in one direction. However if the first 
non-zero a-derivative of V is V"" , then the graph of V may have local minimum or 
maximum in accordance with the sign of V"" . 

(c) Neutral: It may also happen that the potential is identically zero, V(a) = 0, in 
the neighborhood of ao- In this particular case the equilibrium is called neutral 
(throughout this neighborhood). 

In determining the value of V"(ao) recall first that c 2 can be given as 
d0 _ _dQ/da 

Cs " dE ~ dS/da ' 1 ' 

The terms d9/da and dE/da can be expressed by using equations (12.21) and (12.31) . For 
instance, in virtue of (12.5j) . and since a 2 = — 2V(a) and a = — V'(a), we find that dE/da 
can be given in terms of V(a) and V'(a) exclusively. Similarly, dQ/da can be seen to 
depend on V(a), V'(a) and V"(a) exclusively, i.e. no higher order a-derivatives of V(a) 
are needed in evaluating dQ/da. Therefore, in the case of a gravastar in equilibrium, 
i.e. whenever V(a ) = and V'(a ) = 0, V"(ao) can be expressed as a function of the 
variables 77, ao and c 2 = c 2 \ ao . This means — as was mentioned before — in carrying out 
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Figure 1. Typical potentials are indicated. Graph 1 depicts a potential relevant 
for stable equilibrium. The other three graphs, with labels 2, 3 and 4, represent the 
change of the potential caused by dropping more and more massive spherical dust 
shells onto the same initially stable gravastar. In particular, graphs 2 and 3 comes 
with a bounded motion of the resulted gravastar, while in the case of graph 3 the 
exactly maximal mass was dropped onto the initially stable configuration. Graph 4 
corresponds to the situation when the gravastar collapses to a BH. 



the stability analysis of gravastars in equilibrium only the derivative dO/dS = —c 2 of 
the EOS at S comes into play. 

It can be shown that for stable configurations the square of the speed of sound is 
bounded from below by 



2 <99 / da 

'-'s.min 



dH/da 



dQo/dao (2.20) 



on 



where c 2 s min is nothing but the square of the speed of sound c 2 Q at «o for configurations 
with V"(ao) = 0. This, in particular, means that V"{cx.q) > whenever c 2 > c 2 min . 

To see that O2.20p holds note first that S and — given originally as functions 
E(a, d;r/) and Q(a, a, a; rf) — in virtue of the relations a 2 = —2V(a;r]) and a = 
—V'(a;r)) can be given as the functions of the form S = £(«, V(a; 7]); rj) and = 
0(a, V(a; T]), V'(a; rj); rj). Then (I2.19p . along with the relations dE/da = dY,/da + 
(dZ/dV) V and d6/da = dQ/da + (dQ/dV) V + (dQ/dV) V", and the vanishing of 
V'(a;rj) and V"(a;r]) in the case of a gravastar in equilibrium justify the first equality 
in 02 .2 Op . The second one holds because, in virtue of 02. 2p and 02.16P (dT,/da)\ ao = 
(dY,o/da>o), and similarly, in virtue of 02.3p and 02.17P (<96 / da)\ ao = (d®o/da ). By 
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taking into account (I2.16P and ( I2.17P the minimal value of c 2 Q can be given as 

(3-6a +4ag)(l 



4a 2 (l 



-U3/2 



ml? 12 



4(3 - 5a + 2ag)(l - mlf /2 + 8a 2 (a 2 oV - 1)(1 - a l fl 2 ' 



(2.21) 



In figure [2] the configuration space of gravastars in equilibrium is illustrated with the 
help of the dimensionless parameters rj and ceo- In explaining the meaning of some of the 
characteristic curves plotted in figure [2] note that the thin dashed, dashed-dotted and 



;«. 
i'.V. 



0.8 



0.6 



0.4 



lower boundary of de Sitter horizon formation 
upper boundary of positive surface density (E =0) 
upper boundary of DEC (w =1) 
lower boundary of stability region with c s0 =1 
lower boundary of stability region with c s0 =0 
lower boundary of stability region with c s0 =w 




Figure 2. The configuration space of gravastars in equilibrium. The points below 
the uppermost dotted line represent configurations where no cosmological horizon may 
occur in the de Sitter region. At the points of the thick dotted line the surface energy 
density vanishes. Below that line the proper mass of the shell is always positive. 
The thick continuous line is determined by equation (|2.18|) and it shows the upper 
boundary of the region where DEC holds. This curve meets the horizontal axis 
at ao — 25/24 and it has its global maximum at the point M with coordinates 
(ao,w) ~ (1.1503,0.1944). The dashed thick line represents configurations with 
V"(o£q) = and with c^ = 1. Stability holds, i.e. V" > with cj = 1, to the right of 
this curve, while instable configurations are to the left of this curve. The thin dashed 
line represents configurations with V"(cto) = and Cg = 0. This curve also separates 
stable and unstable configurations when c s o is restricted to zero. The dashed thin line 
enters to the DEC region at the point N with coordinates (ao,r)) ps (1.6470,0.0855). 
On the thin dashed-dot line V" = and c^ — — 6o/So holds. The dashed thick line, 
the thin dashed-dot line and the curve determining the boundary of the DEC region 
all intersect at M. Two subregions are indicated by gray shading. In the shaded region 
stable gravastars can be found within the domain where DEC holds and the speed of 
sound at ao does not exceed the speed of light. The dark gray region represents stable 
gravastar configurations irrespective of the value of c!? until it is non-negative. 
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thick dashed lines do correspond to configurations for which the square of the speed of 
sound Cg at a takes the minimal values c 2 min = 0, Wq, 1, respectively. Correspondingly, 
the thick dashed line represents the lower bound of the radius for stable configurations 
for which the sound speed does not exceed the speed of light. If the DEC is also 
required to hold, then ao has to be larger than « 1.15028 which is the ao value of the 
configuration represented by the point M in figure [H 




Figure 3. The configuration space of gravastars on the (So,— @o) plane. An EOS 
may be represented by a curve on this plane. In order to provide a better orientation 
the y = const curves for values ranging from y = to y = 0.18, with separation 
Ay = 0.02, are indicated. The styles of the sheading of the distinguished subregions 
and that of the characteristic curves and points are the same as in figure [T] 

In figure [3] a different representation of the very same configuration space, as 
depicted in figure El is given on the (E , — 0o) plane. The points representing 
gravastars in equilibrium by the values (77, a ) in figure [2] are mapped as (n, a ) 
(Sofa > a o), ©o(?7) «o))- The styles of the sheading of the distinguished subregions and 
those of the characteristic curves and points are in correspondence in these two figures. 
An interesting feature of figure [3] is that, in virtue of the above argument in connection 
with c 2 min , the slope of the tangent at a point of an 77 = const curve is exactly the value 
of c 2 smln relevant for the configuration represented by the pertinent point. 
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3. Dust shell falling onto the gravastar 

The interaction of a gravastar with a spherical dust shell falling onto it will be studied in 
this section. In advance to the arrival of the dust shell onto the surface of the gravastar 
the history of the shell is completely determined by its relativistic equation of motion 
on a Schwarzschild background. To make our investigations somewhat realistic we shall 
assume that the collision of the dust shell and the gravastar is totally inelastic, i.e. all 
the dust particles will be assumed to move together with the surface layer of the yielded 
dynamical gravastar. In determining its basic parameters — the velocity, mass density, 
etc. — the conservation of the 4-momenta will be applied. Finally, by inspecting the new 
potential relevant for the composed system the dynamics of the yielded gravastar can 
be determined. 



d ust sh ell 




M 3 

Figure 4. Left: the de Sitter interior with total mass Mi inside the radius do, 
Schwarzschild exterior with mass parameter M2 gravastar shell with rest mass mi and 
the moving dust shell with rest mass to 2 are shown. Right: the moment of the inelastic 
collision is indicated. 

According to the above-outlined program we shall start with a gravastar in 
equilibrium the total mass of which M2 will be used as our reference unit. The 
dimensionless radius and the mass ratio parameter of this gravastar are = ao/(2M 2 ) 
and 7] = 8/cMf , respectively. We assume that a dust shell, with gravitational mass 
eM 2 , falls onto the surface of the gravastar and the collision between the dust shell and 
the surface of the gravastar is totally inelastic. This, in particular, means that once 
the collision occurs a new three-layer dynamical gravastar forms. Figure H] provides a 
simple depiction of the underlying process. (For another schematic spacetime diagram 
see figure fTUl) Denote the proper mass of the initial shell of the gravastar, of the dust 
shell and of the shell of the gravastar after the collision by m 1; m,2 and m^, respectively. 
The total mass of the system is M3 = (1 + e)M.2. To guarantee to be greater than 
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the Schwarzschild radius of the event horizon of the yielded gravastar, we require the 
inequality e < a — 1 to hold. 

The dynamical characterization of the dust shell in advance to the collision can 
be given based on the results outlined in the appendix. Concerning the velocity of the 
particles of the dust shell at the moment of the collision it will be assumed that either 
it is equal to the velocity that would be acquired by the shell if it was starting to move 
toward the gravastar from rest at infinity, i.e. v 2(00) = 0, or, in the second case, the 
particles of the dust shell will be assumed to be simply attached to the surface of the 
gravastar with zero velocity, i.e. v 2(0:0) = 0@ 

Note that in the 1*2(00) = case the gravitational mass of the dust shell equals to 
its rest mass, i.e. m,2 = sM 2 (see the appendix for more details), and also that, in virtue 
of ( 1A.5I) . the velocity v 2 = dr/dr2 at the moment of the collision can be given as 



v 2 = . (3.22) 

4o 

3.1. Collision 

The conservation of the 4-energy-momentum can be used to determine the radial energy 
and the radial momenta of the surface of the gravastar yielded by the collision. The 
relevant equations are derived in the appendix where a detailed and generic investigation 
of an inelastic collision of two concentric spherically symmetric matter shells is also 
carried out. Because of spherical symmetry the radial component of the 4- momentum 
conservation does correspond to the relativistic 3-momentum conservation — given by 
equation ( 1A.17I) . In virtue of this relation the radial velocity, V3 = da/dr3, of the 
boundary of the yielded gravastar if 1*2(0:0) 7^ can be given as 
87ra ^,o ,£) 

v 3 - v , {6.26) 

^3 

where 



F( s gofa Qq) 2 (1 - 1/oq) (f/2) VlGao + Seap + g 

1 -2 V al + r7 2 ag- 4096 ?r 4 ag S (r?, o ) 4 ' 1 ' 
Here £o(r/, a o) denotes the energy density of the original gravastar in equilibrium given 
by d2HSD. 

The temporal component of the 4-momentum conservation corresponds to the 
energy conservation which can be seen to be equivalent to the dynamical master equation 
of gravastars (12. 2p . Using the parameters of the newly formed gravastar the master 
equation can be given in terms of the above-introduced dimensionless parameters as 



rjal + vj - yjl - (1 + e)/q + vj 



87ra 

6 Although both of these assumptions correspond to certain limiting cases we expect that they provide 
some clue about more generic cases with an initial velocity between two of them. 
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The substitution of (I3.23P into ( I3.25P then yields the following equation: 
£ 3 



\ 



A(r],a ) 2 + 



F(rj,a ,e) 2 



y2 
^3 



\ 



B(a ,e) 2 + 



F(r),a ,e) 2 



Y 2 
^3 



with 



A(r),Ol ) 



l-^o , D , v \A - (1 +e)/a 
— , and //(n n . ej = 



(3.26) 



(3.27) 



£3(77, ao, e) = V A 2 + B 2 - 2VA 2 B 2 + F 2 , (3.28) 

where S3 denotes the energy density of the boundary of the newly formed gravastar. 
Note that the value of E 3 is independent from the EOS. 

Whenever v 2(0:0) = the value of E 3 can be determined by (13.251) where the 
substitution v 3 = has to be used. 



87r«o 57To;o 
Equation (I3.26f) is an algebraic relation for E3 possessing the only positive root 



3.2. After the collision 

As discussed above once the collision occurs E 3 and v 3 can be determined. Note that 
to be able to determine what happens after the collision we also need to know the EOS 
of the newly formed surface. Without this information we cannot determine the new 
surface energy density as a function of the radius, which is essential in assigning the 
new effective potential. 

Concerning the new EOS, since the pressure of the falling dust shell is zero, it 
is tempting to assume that the surface tension remains intact during the process of 
collision. This idea could be supported by a simple-minded summation of partial 
pressures. However, if this happened the particles of the dust should have stopped at 
the surface of the gravastar without interacting with the particles forming the surface 
of the gravastar. This appears to be baseless, and more importantly, to be incompatible 
with the physical picture of inelastic collisions. To resolve the associated discrepancy, 
hereafter we shall assume that instead of the surface tension the EOS of the boundary 
remains intact during the inelastic collision. This means that the dust shell 'melts' into 
the matter of the surface of the gravastar. 

In determining the new potential, as previously, (I2.5P may be applied since after 
the collision a modified three-layer gravastar is formed. The only distinction is that the 
new dynamical variables — indicated by tilde — have to be used in ( 12. 5ft . For the new 
gravastar fj = (1 + e) 2 r] holds as the mass of the outer Schwarzschild region, used as 
a reference unit, changes from M2 to M3 = (1 + e)M2- Accordingly, the dimensionless 
surface density and the dimensionless radius change asa = a/(l + e) and E = (1 + e) E. 
Combining all these observations the new potential may be given as 

V(a; 77) = V((l + e) 2 V , (1 + £)E new (a), a/(l + e)). (3.29) 

where E new (a) is the solution of (12. 8p with the initial condition E new (a ) = E 3 provided 
that the EOS remains intact as assumed. 
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4. Results 

By inspection of the new effective potential the dynamics of the gravastar can be 
determined. If the new potential possesses the form of graph 2 of figure [1] the motion 
of the gravastar surface remains bounded and we can determine the minimum and 
maximum value of the radius. In particular, the maximum mass which may be dropped 
onto the original gravastar without converting it into a BH can also be determined. 

4-1. The applied method 

First note that after collision the initial velocity points inward (unless ^(cto) — 0); 
therefore, the investigation has to start by deriving whether the new effective potential 
has any root in the interval a G [1 , «o] • If not then an event horizon will form and the 
gravastar collapses to a BH. The examination of the potential is made by a numerical 
maximum finding algorithm. If the maximum is found to be less than zero, the collapse 
to a BH will occur. Using an interval bisection method, changing the initial gravitational 
mass of the dust shell to be dropped, the numerical value of the maximal gravitational 
mass, £ max ' M 2 , which may fall onto the gravastar without converting it into a BH, can 
be determined. Whenever W2(«o) 7^ and there is a root in the interval a G [1, ao], the 
motion is bounded and it happens so that a,\ < a < a>2, where a>\ is the largest root 
that is smaller than ao and a 2 is the smallest root that is larger than ao- If ^(cto) = 
then ao coincides with either the upper or lower bound of the motion. 

It is an inherent property of the gravastar model that we cannot have information 
about the de Sitter energy density inside the gravastar, 77, nor about the EOS of the 
shell before and after the collision. To have some quantitative results suitable hypotheses 
have to be applied. In order to avoid the use of a flood of parameters — each of which 
has to possess definite values in numerical investigations — , as indicated in subsection 
12. 3^ in this paper we analyze three types of EOS: homogeneous linear, broken linear and 
polytrop. 

If we fix an EOS with all of its parameters, then, as a function of rj, the stable 
equilibrium radius, a , is to be determined. Note, that this radius can be determined 
by the intersection of the specific EOS curve and the corresponding 7]-level curve in 
the figure [31 Next we have to make an assumption about the EOS after the collision. 
Since we have no observational clues, as it was indicated above, we shall assume that 
the EOS remains intact, which means that the dust shell melts into the matter of the 
gravastar shell yielding an increase of the pressure according to its mass contribution. 
Note that in the case of an initially homogeneous linear EOS the possibility of keeping 
the surface tension to be intact will also be investigated. This latter case might be 
a good approximation whenever the contribution of the dust shell to the new shell's 
composition is negligible. 
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4-2. Homogeneous linear EOS 

As was mentioned in subsection 12. 3|. the simplest possible EOS is the homogeneous 
linear one. If the value of n and Wq are given that of is determined by the implicit 
relation w = — ©o(?7, a )/S (^, a ). It is clearly visible in figure |3] that for a given 
value of wo (= c 2 s £ [0, 1]) the wq So + 0o = line intersects any r\ constant curve at 
most at two points. One of these intersections corresponds to a stable configuration, 
with the smaller So value, while the other is unstable. Accordingly, once the choice 
for i] and Wq is made first we determine the value of «o by numerically solving the 
implicit relation w = — Qq(t}, a )/Ti (r], a ). Then, for any choice of the pair (rj, w ), we 
determine the stable gravastar configuration, and as it was outlined above, the maximal 
mass of the shell that can be dropped onto it without converting it into a BH can also 
be determined. The results of the corresponding investigations are depicted by figure El 
It is worthy to mentioned that some of the special cases depicted on these panels play 
distinguished role in various arguments. For instance, the case wq = 1/3 corresponds 
to ultrarelativistic gas, which — based on certain dimension reduction arguments that 
might not be straightforward to be applied in the present case — is represented in [21] 
by a system with wo = 1/2. 




Figure 5. The normalized value e max / — 1) is shown as a function of r\ and wo for 
the homogeneous linear EOS, where ao stands for the maximal dimensionless radius 
allowed by the DEC. On the top two panels ^2(00) = 0, while on the bottom two ones 
V2(oto) = 0. The panels to the left depict the intact EOS cases, while in the right ones 
the surface tension was kept to be intact. 
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Note that to have a stable configuration the value of Wo = c 2 s has to be larger than 
c smm gi ven by (I2.2ip . In particular, w = c 2 s may take all the values from the interval 
[0, 1] only for i] = 0, i.e. whenever the inner region is vacuum. Also note that e max has 
to be smaller than a — 1. If e > a — 1 occurred, the system composed of the gravastar 
and the dust shell would get beyond its own Schwarzschild radius before the collision 
occurred. This makes the use of the rescalled quantity £ max / (c^o — 1) to be advantageous. 

Before interpreting figures [5] and [6] let us recall that a homogeneous linear EOS, 
6o = —wo S , has wo as the only parameter. Therefore, it is also straightforward to 
investigate the case when both the surface tension and the functional form of the EOS 
remain intact during the collision. It is done by determining the value of the new wq via 
the relation w = — O /S 3 . If this happens the particles of the dust shell do not adopt 
the properties of the particles comprising the surface of the gravastar as only the rest 
mass of the surface will be increased by the collision. In figures [5] and [6] four different 
subcases are considered. These subcases are yielded, on the one hand, by keeping the 
EOS or the surface tension to be intact and, on the other hand, by choosing the velocity 
of the particles of the dust shell at the moment of the collision to be such that either 
1*2(00) = or ^2(0:0) = 0, as discussed in section [3j 
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Figure 6. The wo = 1 sections of the four panels of figure [5] are shown providing a 
better comparison of £ max for the considered configuration. The order of the curves 
from top to bottom coincides with that of the legend at the upper-right corner. 

Figure [6] depict only the wo = 1 — referred to as stiff matter in [8] — sections of the 
four cases indicated in figure [5] which offers a better comparison of the corresponding four 
different dynamical subcases. Figure [7] does also refer to these four basic configurations 
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with Wo — 1. The horizontal sections, corresponding to an 77 = const value, of the grey 
regions indicate the dynamical range of the radius of gravastars. For the particular value 
77 = 0.05 the pertinent new potentials are also indicated. 



0.2 




1 1.5 2 2.5 3 3.5 a 1.5 2 2.5 3 3.5 4 



Figure 7. The dynamical range of the radius for maximally loaded gravastars is 
shown referring to the four basic configurations with wq — 1 considered in figure [6] 
The rj — const sections of the gray regions represent the dynamical range of the radius 
relevant for the considered four subcases. The rj dependence of ao of the gravastars in 
advance to the collision is indicated by the dashed curves. Whenever the initial velocity 
of the dust shell is nonzero (e.g., ^2(00) = 0) the dashed line is in the interior of the 
grey region, whereas whenever V2(cto) = — see, e.g., the right panels — the gravastar 
stars its motion at the edge (minimal or maximal a value) of the dynamical range. For 
the particular choice ij = 0.05 the corresponding potentials, are also plotted, where 
the zeros of the potentials are indicated by the dotted horizontal lines and the applied 
scales are the same on each of these panels. The dotted vertical line with the larger 
a — 9/8 value correspond to the Buchdahl limit, while the other one with a « 1.0384 
indicates the corresponding lower bound to the size of most compact conventional stars 
satisfying the DEC. 

There are two interesting points to be mentioned that are indicated in figure [7J 
First, there are two dotted vertical lines on each of the panels. The one with the larger 
a value corresponds to the Buchdahl limit with a = 9/8 [25], while the other with a 
smaller a value indicates the corresponding lower bound, with a ~ 1.0384, to the size 
of most compact conventional stars. On these panels of figure [7] it is also visible that 
dynamical gravastars cannot get below the Buchdahl limit so they definitely remains 
less compact than the most compact conventional stars satisfying the DEC. 
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Second, on the right two panels with v 2 (ct ) — the grey regions on one side 
are bounded by the dashed lines representing the rj dependence of a of the initial 
gravastar. This is, on one hand, in accordance with intuition as the particles of the dust 
shell are simply simultaneously placed on the surface of the gravastar. What might be 
unexpected, on the other hand, is the following. Whenever the EOS remains intact (see 
the right-top panel of figure [7J there is a subregion where the dashed line bounds from 
above and another where it bounds from belowEI These subregions are separated by a 
point (77, ato) ~ (0.124,1.44) where the initial gravastar remains apparently at rest for 
£ = £ max , while whenever the surface tension remains intact — as it is indicated by the 
bottom-right panel of figure [7J — the dynamical range for e = e max is always bounded 
from above by the initial gravastar configurations. 

4-3. Broken linear EOS 

Since the gravastar radius has always to decrease at the moment of the collision the 
surface density must increase at this moment. Thus, the functional form of the new 
effective potential depends only on that of the EOS relevant for mass density greater 
than S . Motivated by this observation we have used the following particular broken 
linear EOS. So(??,ao) was chosen to be the breaking point, i.e. the choice Si = So was 
made, and only the slope was varied for S > So- Since we have found that the higher 
the value of Wq the more stable is a gravastar, in the broken linear EOS we have chosen 
wi = 1 and W2 to be varied form zero to 1, in steps 0.1, as it is indicated on the right 
panel of figure [8j 

In figure El the r\ dependence of £ max , i.e the maximal value of the gravitational 
mass of the shell that may be dropped onto a stable gravastar without converting it 
into a BH is shown for configurations with 1*2(00) = and ^2(^0) — 0. On both of the 
left panels the value of w 2 decreases from 1 to zero while moving from the top to the 
bottom. It is also visible that for each particular value of w 2 there is a maximal value of 
7] that is indicated by the vertical ticks on the horizontal axis of the the top-left panel. 
These figures do also justify — for the pertinent linear approximation of the EOS — that 
for given rj and ao the larger the value of c s o = w 2 the more stable the gravastar becomes 
in the sense that the allowed value of e max increases together with c s q. 

4.4. Polytrop EOS 

A polytrop EOS of the form of (12.141) has two parameters, A and k. Given a point 
(S , Oo) — with < I ©0 1 < — h is not obvious to choose at all the parameters A and 
k to get an EOS fitting to this point except for the choice k = 2. To overcome this 
slight technical difficulty and also because for a given value of A the choice k = 2 always 

7 The response of the initial gravastar for the placement of the dust shell onto its surface is as follows: 
By increasing the value of e from zero to e max first the radius increases. However, for configurations 
with ao < 1-44 there is a turning point in a and before reaching e max the radius is smaller than the 
initial ag. 
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Figure 8. The applied broken linear EOS is shown on the top-right panel, ao is 
chosen to be maximal such that the DEC is still guaranteed to hold. On the graphs 
of the left panels the value of W2 decreases form 1 to zero, in steps 0.1, while moving 
downward. 

yields the most stable configuration, we have chosen k = 2 and varied only the value of 
A. It is also informative to compare the behavior of a gravastar with a polytrop EOS 
and with a linear one having the same slope at So as it is indicated by the nested panel 
in figure [9] relevant for the particular choice k — 2. For any choice of r\ the value of is 
chosen to be determined by the intersection of the rj = const line and the dashed-dotted 
curve in figure El Note that if «o was chosen to be the maximal value for a given r\ such 
that the DEC holds, the polytrop EOS would degenerate into the homogeneous linear 
one with wq = 1. 

The rj dependence of the maximum of the normalized gravitational mass of the 
dust shell, £ max /(cto — 1)? — where e max is the maximal mass that may be dropped 
onto a gravastar without converting it into a BH — is shown in figure |9] for two types 
of EOS and for the two types of choices concerning the initial velocities, 1*2(00) = 
and V2(qo) = 0. As it might have been anticipated, figure M justifies that a gravastar 
characterized by a polytrop EOS is more stable than the corresponding one with the 
linear EOS in consequence of the fact that the slope of the polytrop EOS is larger than 
that of the linear one for S > S . 
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Figure 9. Right: the applied polytrop EOS with k = 2 and the linear one, having 
the same slope at So a re illustrated. Left: for any choice of r\ the corresponding value 
of ckq is determined by the intersection of the y — const line and the dashed-dotted 
curve in figure [2j The rj dependence of £ max /(ao — 1) is shown for two types of EOS 
and for the two types of initial velocities, ^2(00) = and V2(olq) = 0. The order of the 
curves from top to bottom coincides with that of the legend at the upper right corner. 



5. Conclusions 

In this paper, besides providing a short survey of gravastars in equilibrium, the radial 
stability of gravastars was investigated. The latter was done by determining the response 
of a gravastar to the arrival of a dust shell onto its surface. Concerning the matter 
model of the surface of the gravastar three different types of EOS were applied. We also 
assumed that the DEC holds and the speed of sound c s satisfies the relations < c s < 1. 

Our most important findings are as follows. Among the investigated EOSs the 
homogeneous linear one with c s — 1 appears to provide the largest possible stability 
for gravastars in the sense that for a chosen r\ the value of e max attains its largest 
possible value for c s — 1. It would be interesting to find an analytical justification of 
this observation. It also follows from our investigations that once an EOS is fixed, for 
a given value of oo> the smaller the value of rj, or alternatively, for a given value of r), 
the larger the value of 00 the more stable is a gravastar. The maximal mass of the 
shell that may be dropped onto a gravastar in equilibrium without converting it into a 
BH is also determined. The normalized value of this gravitational mass £ m ax/(oo — 1) 
was found to be the largest e mSiX /(a — 1) ~ 0.96 whenever the matter comprising the 
surface of the gravastar possesses a homogeneous linear EOS with w = 1 and v 2(0:0) = 0. 
Correspondingly the value of £ max ~ 0.96 • (a — 2M)/2M is of the order 1 for a gravastar 
with radius a pa 3M which indicates that gravastars may be significantly more stable 
(although not smaller) than ordinary stars. 

In the appendix a generic framework describing the totally inelastic collision of two 
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shells is also shown. The reason behind considering inelastic collisions is that the de 
Sitter region of a gravastar is unstable against hydrodynamical excitations. Thereby, 
it is reasonable to assume that the particles of the dust shell merge to the particles of 
the surface of the gravastar. In most of the considered cases, for the sake of simplicity, 
we also assume that the EOS of the surface of the gravastar remains intact during the 
collision. All these assumptions indicate that there might be various generalizations of 
the simple model we applied. Nevertheless, our findings in all the investigated cases 
clearly manifest that even a dynamical gravastar cannot be more compact than the 
smallest possible ordinary stars of the same mass provided that for both types of models 
the DEC is guaranteed to hold. This means that the main issue is not that whether 
gravastars can be distinguished from a BH. Assuming that they exist in reality it is 
more appropriate to ask whether they can be distinguished from compact regular stars. 

Let us finally turn the case of a generic BH mimicking gravastar without restricting 
our considerations by referring to the DEC or to any particular EOS of the surface 
layer. Note that in this case the size of the gravastar is not bounded from below so, 
in particular, the assumption a ~ 1 corresponding to the requirement of Mazur and 
Mottola may also hold. Now assume that the mass of this gravastar is M and that 
the gravitational mass of the dust shell — falling onto it — is e ■ M. For simplicity we 
assume that the particles of the dust shell move toward the gravastar with velocity as 
if they started from rest at infinity, i.e. ^2(00) = 0. Then, the mass of the composed 
system is (1 + e)M. In order to avoid the situation in which the composed system gets 
beyond its own Schwarzschild radius Rs = 2(1 + e)M, in advance to the collision, for 
the radius of the initial gravastar, R ao = 2a M, the inequality Rs < R ao must hold. 
This immediately implies then that e < cto — 1, and in turn that the closer the radius 
of the initial gravastar is to the Schwarzschild radius the smaller the mass of the dust 
shell, which may be dropped onto it without converting it into a BH, must be. In other 
words, the better is a gravastar in mimicking a BH the easier is to destroy it by the 
described process. Note that since this argument does not refer to the internal structure 
of the gravastar model it applies not only to the simplified three-layer model of Visser 
and Wiltshire or to the original model of Mazur and Mottola but essentially to any 
analogous construction. 
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This appendix is to introduce the basic equations governing the dynamics of a totally 
inelastic collision of a pair of concentric spherical shells, the outer of which is moving in a 

8 It has been long-known that a gravastar — neither in the shell based model nor in its continuum 
correspondence — can be of considerably smaller size unless the DEC is violated [101 ES] . 
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vacuum spacetime region. These equations are derived by making use of the conservation 
of 4-momentum. In determining them we have applied a suitable adaptation of the 
method of Nakao, Ida and Sugiura which was originally worked out in [27J to describe 
the conservation of 4-momentum in the special case of totally transparent shell crossings. 
In what follows the abstract index notation will be applied. In particular, the lowercase 
Latin indices will indicate the type of tensors while Greek indices will refer to the 
components of tensors with respect to some specified basis fields. 

To start off consider two spherically symmetric concentric shells. Each of these 
shells is assumed to be infinitely thin and their history (before the collision) is assumed 
to be represented by separate timelike hypersurfaces. These hypersurfaces will also be 
referred as shells hereafter. Once the collision occurs it is assumed that the initial shells 
merge such that a single infinitesimally thin shell is formed. As is depicted by the simple 
schematic picture in figure ITUl the spacetime may be divided into three characteristically 
different regions. In advance to the collision shell 1 and shell 2 indicate the interior and 



Shell 3 



Region 1 



Region 3 




Figure 10. A schematic spacetime diagram representing a totally inelastic collision. 
The vertical direction is temporal and time progress upward. 



exterior shells, respectively, while the shell yielded by the merging of these two shells is 
represented by shell 3. The inside spacetime region which has shells 1 and 3 on its the 
boundary is called region 1, the region between shells 1 and 2 is called region 2, while the 
outside region bounded by shells 2 and 3 is called region 3. For notational convenience 
we introduce region 4 which is identical to region 1. Hereafter we use capital Latin 
index K to indicate the label either of a shell or of a region. Accordingly, the value of 
K runs from 1 to 3 unless otherwise stated. 

The metric in the disjoint spherical symmetric regions, labelled by K, will be 
assumed to be given as 

ds 2 K = -f K (r)dt 2 + Mr)-Mr 2 + rW , (A.l) 



Merse Elod Gdspdr and Istvdn Rdcz 



with 



f K (r) = 1 



2M K (r) 



(A.2) 



r 



Note that this class is generic enough to include both the Schwarzschild and de Sitter 
vacuum spacetimes, but the electrovacuum Reisner-Nordstrom solution does also fit to 
this form. 

In restricting our attention to a system consisting of a gravastar and a spherical 
shell falling onto it — the same system investigated in this paper — in regions 2 and 
3, in virtue of Birkhoff's theorem, the metric has to be (locally) isometric to that 
of the Schwarzschild spacetime. Accordingly, in regions 2 and 3 the metric can be 
given by (lA.ljl and (1A.2jl with M^(r) = Mk, where the constants Mk denote the 
corresponding gravitational mass parameters of the pertinent Schwarzschild spacetimes. 
In addition, in region 1 the metric is (locally) isometric to the de Sitter spacetime with 
Mi(r) = (4ir/3)por 3 , where po is the energy density in the de Sitter vacuum. In this 
case, M 3 stands for the total gravitational mass of the system while M 2 denotes the 
total gravitational mass of the entire system. 

Note, however, that in carrying out the following calculations there is no need to 
make any specific choice, i.e. the metric in region 1 will only be required to possess the 
form (1A.1I) and flA.2j) . This provides us the freedom that whenever the interior is chosen 
to be de Sitter type the corresponding equations will immediately be applicable to a 
system formed by a shell falling onto a gravastar, while if the interior is chosen to be 
Schwarzschild type the relevant equations describe the inelastic collision of two confocal 
spherical shells in a Schwarzschild vacuum spacetime. 

In the following denote by u a K the 4-velocity of shell K. The components of the 
4-velocity, with respect to the coordinates in regions bounding shell K, can be given as 



where tk denotes the proper time on the shell K and the suffixes (+) or (— ) indicate 
whether the pertinent region is outside or inside with respect to the shell K. Note that— 
as opposed to the area radius, r, which is everywhere continuous — the time coordinate 
is not continuous across a shell. Once the components of u a K are fixed the spacelike 
outward pointing unit vector, n K , normal to the shell K can be determined by the 
orthogonality condition, n Ka u a K = 0, and its components can be given as 



The evaluation of (I A. 3 1) and (I A. 4 1) requires the determination of the radial velocity, 
drx/drx- There are various ways to get the value of dr K /dr K . The most widely 
known method is based on the thin-shell formalism of Israel-Sen-Lanczos-Darmois. For 
a comprehensive summary of this formalism see [8]. After an appropriate reformulation 




(A.3) 




(A.4) 
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of the 'master equation' (see equation (38) of |S]), which is in fact an energy balance 
equation, dr K /dr K , can be given as 

fdr K V = g K (r) 2 i | M K {r) + M K+1 (r) | m K {rf 
\dTx J mx(r) 2 r 4r 2 ' 

where mjc(f) denotes the rest mass, while gxif) stands for the gravitational mass 
of the shell K. Accordingly, ^(r) = Mjc+xir) — Mk{t) for K = 1,2 and ^(r) = 
M 3 (r) — Mi(r). Note that the rest mass of the shell of a gravastar, in general, is a 
function of the area radius, nevertheless, the rest mass of a dust shell is constant. 

Recall that the components of the 4- velocity cannot be independent as UK a u a K = — 1 
which, in particular, implies that 



dt K{±) = ,J(dr K /dT K ) 2 + f w (r) 

dr K f(±)( r ) 

The substitution of (1A.5j) into (1A.6[) . and some straightforward algebraic manipulations, 
yields that 

dt K (±) _ 9K{r) =f h K (r) 



dr K m K (r)f K{±) (r) ' 
where hxir) = m 2 K (r)/2r, referred as the self-gravity of the shell K. 



(A.7) 



Returning to the main line of the argument recall that our aim is to determine the 
motion of the shell yielded by the merging of the two initial ones. In doing so the only 
guiding principle we can use is the conservation of 4-momentum 

m 3 Uz = m\u\ + m 2 u a 2 , (A. 8) 

which has to be evaluated at the moment of the collision. Note that the content of 
(1A.8I) is nothing but two scalar equations. The temporal component corresponds to the 
energy conservation, while the radial component expresses to the radial momentum 
conservation. The other two components vanish identically due to the spherical 
symmetry. The two non-trivial equations completely determine the motion of the 
outgoing shell since, as we shall see below, both the radial velocity and the rest mass 
777.3 get to be fixed by them. Note that since the collision is inelastic a part of the kinetic 
energy is converted into binding energy which, in particular, implies that m 3 7^ rri\ +m 2 . 

In applying ( 1A.8I) each of the 4-velocities have to be expressed with respect to a 
common coordinate basis field. Start by choosing the coordinate basis field of region 1 as 
our reference frame. Then, in virtue of (IA.3I) . the formal expressions of the components 
of u\ and u% can immediately be determined. Similarly, in virtue of (IA.3I) and (IA.4I) . 
the coordinate basis components of u\ and n\ of region 1 are known, therefore u 2 may 
be expressed there as the linear combination 

14 = -{ulu lb )u$ + {u b 2 n lb )n$. (A.9) 

By applying (jSUD, (D, and JQD 
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Ur(±)c 



dr K d*K-(±) 



,0,0 , 



dr K dr K / 

, / ,dt K (±) dr K /dr K 
-/(±)(" 



0,0 



dr^ /(±)(r) / 
which, along with the notation px = /^(dr^/drR-), implies that 



(A.10) 
(All) 



u 2 u lb 



u 2 n lb 



PlP2 ~ (ffl ~ fol)(fl2 + M 

mim 2 / 2 
(ffi ~ ^jjPg ~ (gg + h 2 )pi 
m 1 m 2 f 2 



(A.12) 
(A.13) 



By substituting these relations into (1A.9j) the components of w 2 in region 1 are 



determined. Finally, by making use of these components, along with (1A.7|) . the two 
non-trivial equations in (IA.8j) can be given as 



h 3 -g 3 = h x -g x - 



(92 + h 2 )(pl - g\ + hi) + 2 Pl h lP2 



P3 = Pi- 



m 2 J 2 

p 2 (pl ~ gf + hp + 2p 1 h 1 (g 2 + h 2 ) 
mlf 2 



(A.14) 
(A.15) 



where g 3 can be eliminated from (lA.14j) by using the relation g 3 = gi + g 2 - We would 
like to emphasize that in deriving these relations no use of the EOS of the shells has 
been made. 

Let us finally restrict our considerations to the case of a dust shell falling onto a 
gravastar. Since the gravastar is in equilibrium pi = 0. Thus, in virtue of (IA.15j) . we 

get 

hi 



m 3 v 3 = m 2 v 2 



9l 



(A.16) 



Whenever v 2 = then v 3 = also holds. Whenever v 2 7^ and the gravitational mass 
of the shell is equal to its proper mass — in this case v 2 (oo) = — by making use of the 
notation of section 3 the relations m 2 = eM 2 , m 3 = 47r(a 2M 2 ) 2 £ 3 /M 2 , f 2 = l — l/a , 
gi = M 2 (l — r/al) and (I3.22p can also hold. By substituting them into (1A.16j) the radial 
velocity, v 3 , of the shell, at the moment of the collision, can be written as 



^3 



4vr £g ( ao _ 1) e y / £ 2 + g £QQ + iq aQ 
£ 3 (1 - 2 77 a[j + rfal - 4096 7r 4 a(j Eg) 



(A.17) 
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